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Abstract 



We give a framework for molecular multiphoton excitation process induced by 
an optimally designed electric field. The molecule is initially prepared in a coherent su- 
perposition state of two of its eigenfunctions. The relative phase of the two superposed 
eigenf unctions has been shown to control the optimally designed electric field which trig- 
gers the multiphoton excitation in the molecule. This brings forth flexibility in designing 
the optimal field in the laboratory by suitably tuning the molecular phase and hence by 
choosing the most favorable interfering routes that the system follows to reach the target. 
We follow the quantum fluid dynamical formulation for designing the electric field with 
application to HBr molecule. 



I. INTRODUCTION 

Finding ways to control the quantum phenomena is, arguably, the essence 
of chemical and physical processes. Advancement in generating tailored light pulses has 
made it possible to steer matter towards a specific goal, namely, to control its future [1-7]. 
It is now well established that light shaping is an efficient way of controlling quantum 
systems and thus, for example, chemical reactions. Recent days have suggested a variety of 
control schemes and their implementations in the laboratory [8-14] based on the coherent 
nature of the laser radiation. 

Coherent control (CC) scheme accesses the broad range of quantum interference 
effect through the relative phase of the light waves and/or the relative phase of the 
molecular eigenfunctions as has been theoretically demonstrated by Brumer and Shapiro 
[5,6] and experimentally realized by others [14]. A somewhat different control scheme 
uses the rapid progress in ultrashort laser pulse technology, e.g. , Tannor et al [2,15] 
suggested a "pump-dump" technique, which has been realized experimentally by several 
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groups [16,17]. These techniques deliver control through a limited number of parameters, 
viz., the optical phase, the molecular phase and the time delay between pump and dump 
laser pulses. 

In the general field of controlling dynamics these parameters are not sufficient and 
it was Rabitz et. al [18] who first proposed in their optimal control theory (OCT) that the 
electric field of the laser pulse should be treated as parameter and be specifically designed 
both temporally and spectrally. Theoretically, OCT creates a nonstationary state of one's 
choice, by optimally designing the control field. The control field couples the system with 
itself in a way that as the system evolves, its motion is modified by the influence of the 
field along its path, and the optimal control field provides the maximum overlap of the 
time-evolved state with the state of one's choice. Although there have been progress in 
the pulse-shaping technology, it still remains a huge task ahead to design such field in the 
laboratory. 

In this paper, we present a theoretical treatment of the OCT which introduces the 
parameter, molecular phase, mostly used in the CC scheme. This prepares the molecule 
in a coherent superposition state comprising of two molecular eigenfunctions, prior to its 
being submitted to the control field. The molecular phase is shown to be an experimental 
parameter which can take a whole range of values for designing the control field. We orga- 
nize this paper as follows : section 2 describes the preparation of the initial superposition 
state comprising of two vibrational eigenfunctions of the ground electronic state of HBr 
molecule. This can be accomplished by non-resonant two-photon interaction between the 
molecule in the ground state and two laser electric fields. The relative phase of the laser 
fields defines the molecular phase of the eigenfunctions forming the superposition. Section 
3 briefly presents OCT within the framework of quantum fluid dynamics (QFD), called 
OCT-QFD method, for desiging an electric field to meet a specific objective. The details 
of the OCT-QFD have been published before [19,20]. Section 4 applies the OCT-QFD 
for the manipulation of the multiphoton excitation of HBr molecule subject to the design 
of an optimal control field. Section 5 concludes the paper. 
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II. MOLECULAR COHERENCE AND SUPERPOSITION STATE 



Here we showed one possible way to prepare the superposition state of HBr. The approach 
is based on a nonresonant two-photon transition phenomena whereby the molecule intially 
in the eigenstate \v\J\M\ > transfers population to the eigenstate |z/ 2 J 2 M 2 > through the 
intermediate states. Here \vJM > refers to the ro- vibrational eigenfunction correspond- 
ing to an electronic state, where v is the vibrational quantum number, J the rotational 
quantum number and M the projection of J onto the internuclear axis in the laboratory 
coordinate. Let us consider the superposition where both \u 1 J 1 M 1 > and |z/ 2 J 2 M 2 > 
belong to the same electronic state, viz., the ground state 1 £ + of HBr molecule. The 
two-photon excitation is assumed take place by the laser field given by 

E(t) = fjE^f(t)cosu L t(l + cos<pl) (1) 

where 

f(t) = e-^> 2 (2) 

with t being the switch-on time of the pulse, r being related to the full width at half 
maximum(FWFM) as FWHM = 2r^J —ln(l/2), uj l the central frequency of the pulse 
and Ef ] the maximum amplitude of the laser electric field. The electric field polarization 
vector rj is considered along z direction in the laboratory coordinate. This field can be 
obtained by the superposition of two identical laser pulses with relative phase 4>l- The 
superposition state is thus given by 

= c 1 (t)\iy 1 J 1 M 1 > e' Wlt + c 2 (t)|z/ 2 J 2 M 2 > e~ W2t (3) 

At t — —00 |cx| 2 = 1 and at time t |ci(t)| 2 = 1 — |c 2 (t)| 2 . In the above equation loi = ^, 
where E i is the energy of the \v i J i M i > eigenstate. For two-photon interaction, applying 
the second order time- dependent perturbation theory we obtain 

c 2 (t) = -t^ci^cM, /* e^-^+^W E{t')G(t')dt' (4) 
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where G(t') is given as 

G{t') = f j^-W'+WW' E(f)df (5) 

I refers to the intermediate state involved in the two-photon absorption process, IYs are 
the life time of the ith state and djj is the dipole matrix elements between i and j states. 
The summation involving only the dipole matrix elements in Eq.(4) can be written as [21] 

E M = £(-i) M+Ml -*- Kl ^ (6) 

F ! , m (2J 1 + l) 1 / 2 (2J 2 + l) 1/2 4/ ) (^) 



Ji J -h 
-Mi - M M-> 



1 Jl J J 2 ^ 

v -K x -K K 2 j 




The function F UlUl is the Frank Condon factor whereas the term (rj) is purely geometric 
depending on the direction of the polarization vector of the laser field and is given by 

L^( V ) = (2J+lf^(-^) A+B VAVB (7) 
/ 

V 

where A and B run over -1, and 1 corresponding to X, Z and Y respectively in the 
laboratory coordinate. The term M^\ii e gI ) is free from the experimental conditions which 
is purely molecular and is given by 

Mi J) (/4,) = (2</+l) 1/2 E(-l) a+6 (8) 



where a and b run over -1, and 1 corresponding to x, z and y respectively in the 
molecular coordinate. The function (n e gI ) a = (A*f fl )a is the a-th component of the elctronic 
dipole matrix element between the ground electronic state and the intermediate states. 
It is obvious that the intermediate state must be either S + (when the transition dipole 
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operator is — J2i z i [22]) or II (when the dipole operator is ^= J2i(%i±Wi) [22]). Whereas Ji, 
J 2 , J, K u K 2 , K, Mi, M 2 and M take the values as K=0, M=0, M 2 = M 1 , K X = K 2 = 0, 
J 2 = J\ (or J 2 = Jl + 2). The vibrational quantum numbers for the superposition can be 
either v 2 = v\ or u 2 ^ V\. Note that the assumption of the perturbation theory implies 
that |c 2 (t)| 2 cannot exceed = 0.2 which, in other words, restricts the power of the pump 
laser [23]. Notice that the molecular phase defined as <\>m = ^ an_1 (§(^|y> where 3? and S 
refer to the real and the imaginary parts respectively, depends mainly on <f> L , lo\ and u> 2 . 
The factor (Tj + T 2 )/2 or (Tj + Fi)/2 in the exponent of equations 4 and 5 respectively, 
have the least contribution to the molecular phase since they are neglible compared to 
uj 2 i = uj 2 — uji or uji = uji — uj\ respectively. Thus one can vary 0m in the laboratory by 
varying <f>L, Ui and u> 2 . 



III. OCT-QFD METHOD 

In optimal control theory an objective functional, corresponding to a specific dynamics of 
one's choice, is minimized with respect to the electric field by solving the time-dependent 
Schroedinger (TDSE) equation. Consider a general target expectation value defined as 

6 T = [ T Qp(x,T)dx (9) 

Jo 

where 6 is an observable operator and p(x,T) = ty*(x,T)ty(x,T) with ^(x, t) being the 
complex wave function at the target time t=T. This wave function ty(x,t) obeys the 
TDSE 

^-^f 1 = h^V 2 + V(x) + V ext (x, t)Mx, t) (10) 

where V typically confines the particle in a locale and V ext is the control taken here as 
—fi(x)E c (t) with E c (t) being the control electric field to be designed and /jl(x) the dipole 
moment. 

The goal is to steer Ot as close as possible to a desired value Q d . The active spatial 
control interval is taken as x\ < x < x r over the time < t < T that the control process 
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occurs. We now desire to minimize the cost functional J cost = Jtarget + J field, where Jtarget 
and J fi e i d are given by 

1 



Jtarget = ^ X (Q T - Q d Y (11) 

and 

J field =\"e [ E 2 c (t)dt (12) 

The minimization of J cost with respect to E c (t) must be subjected to the satisfaction of the 
equations of motion for ^(x, t) in Eq.(lO), which can be transformed into two equations, 
viz., the continuity equation 

^ + V-(pv) = (13) 

and a modified Hamilton Jacobi equation 

*> + ™^ + V + V ^ + V. = (U) 
at 2m 

with the substituition ^(x,t) = p 1 / 2 (x,t)e iS{x ' t)/h in Eq.(lO) where V q = ~^~4rr- = 
-^[V^np 1 / 2 + (V/np 1 / 2 ) 2 ]. This forms the basis of the QFD [24,25] treatment of TDSE. 
This equation can be transformed into the one for the evolution of the velocity vector v 
by taking the gradient to give 

A v = -(v • V)v - -V(V + V ext + V q ) (15) 

Defining the quantum current as 

j(x, t) = --9f[**(x, t) V*(x, *)] = p(x, t)v(x, t), 
m 

one readily obtains the equation of motion for j as 

jrj = -v(V • j) - (j • V)v - ^(W + + V q ) (16) 
ot m 

Thus within the QFD formulation, we need to minimize J cost with respect to E c (t) subject 

to the satisfaction of the equations of motion for p and j given by Eqs.(13) and (16) 

respectively. 

We may fulfill this constraint by introducing the unconstrained cost functional as 
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where Ai(x, t) and A 2 (:r, £) are Lagrange's multiplier functions. 



An optimal solution satisfies 5 J = 0, which is assured by setting each of the functional 
derivatives with respect to Ai, A 2 , p, j and E c to zero. The first two, i.e., the functional 
derivatives with respect to Ai and A 2 regenerate the QFD equations viz., Eqs.(13) and 
(16). The three others are obtained in the forms : 

^ + £(A 2 ^) + 5 1 [p,j,A 2 ] = (18) 
' + tHAiUa) - X ^( V + V q (X 2 ) + V ext ) + 5 2 [p, j, A 2 ] = (19) 



and 



where 



and 



dt dx dx 



5E c (t) 



SJ r Xr d 

/ X 2 (x,t)p(x,t) — p(x)dx + uj e E c (t) = (20) 

J xi OX 



a = 2^ (2i) 



A 2 9 . . . j 2 d\ 



S2 = -~(V q (p) - V q {\ 2 )) - (22) 
m ox p 2 ox 

h 2 d 2 r 1 d 



4 m 2pl/2 Q x 2 pl/2 Q x 

d 2 1/2 d_ 

4 m 2p3/2 Q x 2 Q x 



v ^> = = ~h vHnXr + (v " a > /2)21 (23) 



The corresponding final conditions are 
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u x [Qt - Q d Mx) - X^x^) = 



(24) 



and 

A 2 (a;,T) = (25) 

The equations (18) and (19) for A 2 and Ai respectively ressemble that of p and j with the 
only difference being the extra source terms Si and S 2 . The source terms depend on p 
and j. v\ in the above equations is the 'velocity' associated with the Lagrange's multiplier 
and is given as v \ = There are now two different quantum potential terms, one of 
which is a function of p(x,t) and the other is a function of X 2 (x, t). In this formalism 
the evolution of Xi(x,t) takes place by V q (X 2 ) as well as the difference of the two types 
of quantum potential. In obtaining the above equations we have standardly assumed no 
variation of either p(x, 0) or j(x,0). Thus, we start from the initial value of p(x,0) and 
j(x,0) to solve Eqs.(13) and (16). Eqs.(18) and (19) can be solved for X 2 (x, t) and Xi(x, t) 
by integrating backward from time T using X±(x, T) and X 2 (x,T) given in Eqs.(24) and 
(25) respectively. The equations (13), (16), (18) and (19) are non-linear thereby calling 
for iteration to solve(see ref. [24] and [25] for details). Finally the desired control electric 
field is given from Eq.(20) as 

E c (t) = -— r X 2 (x, t)p(x, t)^-fi(x)dx (26) 

Ul e Jxi OX 

4 APPLICATION TO HBR MOLECULE 

The OCT-QFD has been applied for manipulating the multiphoton excitation process of 
HBr molecule whose initial density p(x, 0) is given by 

p(x,0) = |ci| 2 pi(x, 0) + \c 2 \ 2 p 2 (x,0) +4|ci||c 2 |cos0 M pi2(a;, 0) (27) 

where c\ and c 2 have been obtained following section 2. In Eq.(27) Pi(x,0) = \ipi(x, 0)| 2 
and the pi 2 (x, 0) = ipi(x, 0)ip 2 (x, 0) where ipi(x, 0) corresponds to the \uiJjMi > eigenstate 
with x being the internuclear distant. Although the perturbation theory permits |c 2 | 2 < 
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0.2 there could be other non-perturbative methods resulting |c 2 | 2 > 0.2. Thus, in the 
results below whenever we consider | C2 1 2 > 0.2 we assume the existence of some non- 
perturbative methods. A whole range of <pM can be attempted to follow the different 
interfering routes to reach to the target state by suitably modifying the control field 
E c (t). 

Equation (4) for c 2 shows that we need the ro- vibrational eigenfunctions in the ground and 
the excited electronic states. These have been evaluated by solving the time independent 
Schrodinger equation using the Fourier Grid Hamiltonian(FGH) method [26]. The ground 
and the excited electronic potentail energies have been taken from the ref. [27]. Figure 
1 shows c 2 for the superposition ci|000 > e~ iElt/h + c 2 |100 > e ~ iE2t/h where E 1 and E 2 
are the energies of the eigenstates |000 > and 1 100 > respectively. This gives |c 2 | 2 = 0.2, 
| ci | 2 = 0.8 and 4>m = 0.25 radian after the laser fields disappear. For any superposition 
one can choose a whole range of values for <pM for a given |ci| 2 and |c 2 | 2 by choosing 
different values of 4>l- This (pM, Wi, u 2 , \ci\ and |c 2 | have been used as parameters for 
designing the control field for the occurrence of a specific dynamical change in the molecule 
corresponding to the target operator 9 = x, where x represents the average distant along 
the internclear axis. Other kinds of target operator can also be considered in the present 
QFD-OCT following some modifications of the control equations (see ref. [20]). In the 
results shown below the spatial range of the calculation was < x(a.u.) < 14 and the 
time interval was < t < T with T=5000 a.u. The total number of the spatial mesh 
points taken was 52 which gives 5x = 0.27 a.u. whereas the total number of time steps 
was 2000 which gives St = 2.5 a.u. The weight uo e is taken 0.5 and oo x 500 and the desired 
target value Q d = 3.2 a.u. 

Figures 2,3,4 and 5 show the results for the optimal control field and the dynamics in- 
volved. Notice that the desired target, i.e., Q d = 3.2 a.u. is the same for alll the cases 
(Figs. 3, 4(b) and 5((b),(d))) although the paths through which the dynamics occur to 
reach the desired target were different. Basically, the phenomena that causes the ex- 
pectation value of x to change in time is the multi-photon excitation of the molecule 
induced by the control field along a selective path. If we assume that the control field 
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strengths are low enough for the perturbation theory to be valid, which is indeed true 
in the present results(e.g., E max pa 0.3 a. u. (Fig. 2(c)) corresponding to maximum inten- 
sity ~ 10 14 W att / cm 2 and E min pa 0.06 a. u. (Fig. 5(c)) corresponding to pa 10 13 Watt/cm 2 
intensity) we find that the time evolved wave function under the control field is 



where the summation j occurs over all the possible states excited by the optimal field 
starting from the superposition state. In the above equation Cj(t) depends on C\ and c 2 
as Cj(t) = Cifi(t) + c 2 f 2 (t), e.g., in the first order perturbation theory fi, is given as 



where E(t) is the optimal control field. Thus, the density at any time after the control 
field is on is given by p(x,t) = p ni (x,t) + p in (x,t) where "in" referes to the interference 
term and "ni" to the non-interference term. Similarly the target expectation has both 
non-interference and interference terms in it, i.e., < x > (£) =< x >„« (t)+ < x > in (t). 
This readily gives the optimal control field E(t) = E ni (t) + E in (t) The non-interference 
terms associated with p ni (x,t), < x > ni (t) or E ni (t) have two parts, one relates to |ci| 2 
and the other to | C2 1 2 whereas the interference terms depend on |ci|, | C2 j and 0m- 
Thus, one actually controls the dynamics of the multiphoton excitation process subject to 
the suitable designing of an optimal electric field which itself can be controlled by varying 
the parameters, viz., 0m, |ci|, | C2 1 , 0J\ and u>2- This has been depicted through the figures 
2 to 5. The dynamics by the control field, take place through three routes : (a) route that 
ends up at ipj (Eq.(28)) from \v 1 J 1 M 1 >, the probability of which is proportional to |cx| 2 , 
(b) route that ends up at ipj from \u2J2M2 >, the probability of which is proportional to 
I C2 1 2 and (c) the route that end up at ipj through the interference between the routes (a) 
and (b), the probability of which is proportional to |ci||c 2 |sm0M or l c i| |c2|cos0m- Thus 
the molecular phase can cause certain excitations by the control field during the process 
of its designing, which are not present in the absence of the superposition state. 
Figure 2 shows the control fields corresponding to four different values of M for the 
superposition |000 > +|100 >. These fields excite several vibrational excited states(not 



(28) 



j 




(29) 
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shown here) in the process of achieving the target state, i.e., a state with <d d = 3.0 a.u. 
The peak value of the strongest field (Fig. 2c) is « 0.2 a.u. (corresponding intensity is 
~ 10 14 W att/cm 2 ) which can be readily attained in the laboratory. However, the pulse 
shape presents a chalange to the present-day laser shaping technology. A detailed char- 
acterization of the optimal field can however, be made by Fourier transforming the fields. 
Fig. 3 shows the average distance < x > as a function of time, the corresponding control 
fields are shown in Fig.2. Notice that the molecular phase, 0m changes the course of 
the dynamics of excitation as evidenced by the behaviour of the the expectation value of 
< x > and the corresponding optimal control fields. 

Figure 4 shows the optimal electric field and the corresponding expectation value of x for 
the superposition |000 > +|200 >. This result can be compared to that of Figs. 2a and 
3a for the superposition |000 > +|100 >. Figure 5 on the otherhand, shows < x > (t) 
and the corresponding optimal field for different values of | C2 1 2 for the superposition 
|000 > +|100 >. 

IV. CONCLUSION 

This paper presents the optimal design of the electric field by using the QFD formulation. 
Molecular coherence has been introduced in such design of the electric field, by creating 
a nonstationary superposition state comprising of two vibrational eigenfunctions of HBr 
molecule in its ground electronic configuration, prior to its submission to the electric 
field to be designed. The molecular coherence is created by allowing the molecule to 
interact with a superposed laser field. We showed that the molecular phase <pM can be 
experimentally varied and used as a parameter to modify the optimal electric field so as 
to manipulate certain dynamical change in the molecule. Applications to the multiphoton 
excitation of HBr molecule, described by the expectation value of x, show that control 
range is extensive. The results show that pulses with different structures result in different 
excitation processes since they follow different interfering routes, and the interference 
routes are controlled by the molecular phase. One can indeed optimize other quantity 
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with some constraints in it so as to follow only one out of several interference routes. 
The extent to which molecular coherence (and hence, the quantum interference) enters 
into the optimal results of the elctric field is central to the understanding of the control of 
the excitation. At present we found a variety of results. In several cases such interference 
was indispensable in producing the optimal results, often necessitaing the involvement 
of large numbers of interferring routes to the excited states. Although the multiphoton 
excitation phenomena are not so suitably described by the operator x, this is the first 
case study where we combined the coherent control method with the optimal control 
method. Other multiphoton excitation processes, with suitable operators, are currently 
under investigation. 

Note that one could use the standard control equations based on the TDSE, however, re- 
cent studies [20,24,?,?] show that the QFD equations require lesser number of spatial grid 
points than the TDSE when solved numerically which apparently enhances the efficiency 
and the numerical saving of the QFD-OCT. This enhancement is attributed mainly to 
the relatively slow varying nature of the hydrodynamical variables, compared to the wave 
function, in the spatial coordinate. It may be mentioned that for other target operators 
the TDSE based control equations may be easier to handel than the QFD based equations 
however, in the present study the chosen operator is such that QFD based equations are 
more easy to solve. 
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Figure Captions 

Figure 1 : Time(in picosecond) variation of |c 2 (t)| 2 (a), K(c 2 (t)) (b), ^s(c 2 (t)) (c) and 
0m(£) (d) corresponding to ci|000 > e~ Wl * + c 2 (t)|100 > e~ U2t superposition state, where 
K and S refer to real and imaginary respectively. This superposition is obtained with a 
short (spectrally) laser pulse whose central frequency corresponds to 432.61 nm, FWHM 
corresponds to 5.54 ps = 5.31 cmT 1 t = and the peak intensity is 2.647 x 10 6 Watt/cm 2 . 

Figure 2: Optimal electric pulse, E(t) in atomic unit (a.u.) plotted against time in a.u. 
Label (a), (b), (c) and (d) correspond to the molecular phase , 0m = tt/5, — 7r/5, n/3 and 
7r radian respectively. The superposition state is ci|000 > e~ lElt l h + C2 (t) 1 100 > e~' lE2t ^ h 
where | c x | 2 = 0.8 and |c 2 | 2 = 0.2 and Q d = 3.2 a.u. 

Figure 3: Average distant, < x > in atomic unit (a.u.) plotted against time in a.u. Label 
(a), (b), (c) and (d) correspond to the molecular phase , 0m = 7r/5, -tt/5, 7t/3 and 7r 
radian respectively. The superposition state is ci|000 > e ~ lElt / h + c 2 (t)|100 > e ~ lE2t l n 
where |ci| 2 = 0.8 and |c 2 | 2 = 0.2 and Q d = 3.2 a.u. The corresponding electric fields are 
shown in Fig. 2. 

Figure 4: Optimal electric field (label (a)) and average distant (label (b)) in atomic 
unit (a.u.) plotted against time in a.u. for <p M = ir/5. The superposition state is 
cx|000 > e~ iE ^ h + c 2 (t)\200 > e - iE ^ h where | c x | 2 = 0.8 and |c 2 | 2 = 0.2 and Q d = 3.2 a.u. 

Figure 5: Optimal electric field (label (a) and (c)) and the corresponding average distant 
(label (b) and (d)) plotted against time for 0m = tt/5 radian and <d d = 3.2 a.u. The 
superposition state is Ci|000 > e~ iElt/n + c 2 (t)|100 > e~ iE2t/n where |c 2 | 2 =0.3 (label (a) 
and (b)) and 0.1 (label (c) and (d)). 
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FIG. 2. 
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FIG. 4. 
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